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ABSTRACT 

A moderately massive early Sun has been proposed to resolve the so-called faint 
early Sun paradox. We calculate the time-evolution of the solar mass that would be re- 
quired by this hypothesis, using a simple parametrized energy-balance model for Earth's 
climate. Our calculations show that the solar mass loss rate would need to have been 
2-3 orders of magnitude higher than present for a time on the order of ~ 2 Gy. Such 
a mass loss history is significantly at variance (both in timescale and in the magnitude 
of the mass loss rates) with that inferred from astronomical observations of mass loss 
in younger solar analogues. While suggestive, the astronomical data cannot completely 
rule out the possibility that the Sun had the required mass loss history; therefore, we 
also examine the effects of the hypothetical historical solar mass loss on orbital dynam- 
ics in the solar system, with a view to identifying additional tests of the hypothesis. 
We find that ratios of planetary orbital spacings remain unchanged, relative locations of 
planetary mean motion and secular resonances remain unchanged, but resonance widths 
and the sizes of the Hill spheres of all planets increase as the Sun loses mass. The pop- 
ulations and dynamics of objects near resonances with the planets as well as those of 
distant irregular satellites of the giant planets may contain the signature of a more mas- 
sive early Sun. However, this would not be an unambiguous test, as other dynamical 
processes in the early Solar System, such as the growth and presumed migration of the 
gas giant planets, may also be adequate to explain such signatures. 

Subject headings: Earth — solar system: general — solar-terrestrial relations — stars: 
mass loss — Sun: evolution 



1. Introduction 



It has been well established in models of solar evolution that solar luminosity has increased 
over the lifetime of the Sun due to the increase in the mean density of the solar core as hydrogen 
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is converted into helium ([Goughl Il98ll ). This process has led to a ~30% increase in the solar 
luminosity during the past 4.56 Gy. Because the Sun's luminosity was significantly lower in the 
past, calculations of Earth's ancient climate based on contemporary values of terrestrial albedo 
and atmospheric composition suggest that the mean surfa ce temperature of Earth would have been 
below the freezing point of seawater until roughly 2 Ga (jSagan &; Mullen! 1 19721 ) . However, there 



is ab undant geological evidence that Earth had surface liquid water as early as 3.8 Ga (jKasting 
1989). Geological evidence also suggests that Earth's climate may even have been warmer on 
average than the contemporary climate, with t he exception of the relatively brief snowball Earth 
epochs at ~ 2.2 Ga, ~ 7 60 Ma, and ~ 620 Ma (|Hvde et al-lboool : iHoffman et al1ll998l : iKirschvink 



1992 



Evans et al.lll997l ). The problem is even more pronounced for Mars, where, even at the 



current levels of solar luminosity, the mean surface temperature of Mars is far too cold for liquid 
water to exist, and yet there is abundant evidence that Mars was war m enough for surface liquid - 



water and that surface liquid water was present at around 3.8 Ga (jGoldspiel Squvres 



Jakoskv & Phillipsll200lh . 



1991 



The contradiction between the cold ancient terrestrial and martian climates predicted by the 
lower solar luminosity and the warm ancient terrestrial and martian climates derived from geolog- 
ical evidence has been called the faint early Sun paradox. Most attempts to resolve the paradox 
involve models of the earl y atmospheres of Earth and Mars containing enhancements of atmo- 
spheric greenhouse gases. ISagan fe Mullen! (|1972l ) proposed an early atmosphere of Earth that 
contained small amounts of ammonia (NH3) to provide enough of a greenhouse effect to offset the 
lower solar luminosity. However, it has been shown that NH3 would have been rapidly dissociated 
by sol ar UV radiation, and was unlikely to be a major cons tituent of the early atmosphere of 
Earth (jKuhn fc Atrevalll97fll : fKastinglll982 : IPavlov et al.lboOll ). 



Carbon dioxide (CO2) has also been proposed as a greenhouse gas capable of resolving the 
paradox. A CO2 concentration of at least 100 x present atmospher ic level (PAL) would have to 
have been prese nt to prevent the Archean Earth from freezing over (jKasting! Il987l ). According to 



Rye et al.l (1 19951 ). this abundance of CO2 in the atmosphere should have led to the common presence 
the mineral siderite (FeCOs) in Archean paleosoils, but the absence of siderite from known soils 
older than 2.2 Ga suggests that CO2 could not have been present in the Earth's atmosphere in the 
required amount. Analysis of the weathering rinds of river gravels dated to 3.2 Ga suggest that 
the presence of Fe(II)-rich carbonate in the rinds sets a lower limit of CO2 partial pressure in the 
atmosphere to only several times the present value, which is two orders of magnitude belo w what 
is requir ed to keep the surface temperature of the Earth above freezing (jHessler et al.ll2004l ). More 
recently, lOhmoto et al.l (J2004J) have challenged the interpretation of Rye et al. regarding the lack 
of siderite in paleosoils, arguing that if even a very small amount of oxygen was present in the 
Archean Earth's atmosphere, siderite (or any other ferrous-rich mineral) would have been unstable, 
and CO2 could still have been present in high enough concentration to warm the Earth. 

Methane (CH4) has a lso been proposed as a possible gree nhouse gas responsible for keeping 
the Archean Earth warm (|Hartlll978l ; iKiehl fc Dickinson! Il987l ). Although CH 4 , like NH 4 , is sus- 
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ceptible to photodissociation by solar UV, it can remain in the atmos phere for much lon ger and 
methanogenic bacteria could maintain the required atmospheric levels (jPavlov et al.l l2000l ) . Were 
there adequate populations of methanogenic bacteria available in the early Archean? With only 
trace amounts of evidence of life even existing prior to about 3.6 Ga, it is difficult to support the 
evidence of a warm early Earth throughout the Hadean and early Archean with methanogenic 
bacteria replenishing atmospheric CH4. 

For Mars, one long standing problem with solving the faint young sun paradox with a dense 
CO2 atmosphere is that models predict that the CO2 will begin to condense out as a cloud layer 
at the required atmospher ic pressures, th ereby increasing the global al bedo and off-setting the 
warming greenhouse effect ( Kastinjll99ll ). IPierrehumbert Erlick (1998) proposed a model of the 
early martian atmosphere in which CO2 cloud particles of a certain size could cause a scattering 
greenhouse, requiring an atmospheric pressure of less than 1 bar CO2 to keep early Mars above 
the freezing point of water. However, laboratory studies of CO2 cloud formation under martian 
conditions seem to suggest that the types of clouds that could form on Mars, even with a CO2 
atmosphe re with pressures as high as 5 bars would not warm the planet above the freezing point 



of water (Glandorf et al. 



2002 



Colaprete fc Toonl 12003 1. 



An alternative solution to the faint early Sun paradox involves a non-standard solar model in 
which the Sun has lost significant mass over time. A more massive early Sun will have two effects. 
First, since stellar luminosity is mass-dependent, a larger solar mass implies a correspondingly 
larger solar energy output. Second, owing to the existence of adiabatic invariants of the keplerian 
orbits, the planets would have orbited closer to the Sun had the solar mass been higher. There 
is an upper limit to how much more massive the Sun could have been in the past based on the 
constraint that if the solar flux at Earth had been ~ 10% higher at any time in the past then the 
Earth would have lost its water due to a moist greenhouse atmosphere, in which water re aches the 
stratosphere and is lost due to UV dissociati on and escape of hydrog en (jKasting] Il988l ). Such a 
process is thought to have occurred on Venus (jKasting Pollack! Il983 ). 



Whitmire et al.1 ()1995l ) suggested that a young Sun 3%-7% more massive would be able to 



explai n the evidence of liquid water on the surface of Mars at 3.8 Ga. ISackmann Boothroyd 
(j2003 ) tested solar models in which the Sun first entered the main sequence with an initial mass of 
up to 1.07 Mq. Using various mass loss rate functions that were constrained such that after 4.56 Gy 
on the main sequence the solar models had the currently observed solar wind mass loss rate and 
the current mass for the sun, they compared the resulting solar interiors of their models against 
solar interior profiles (adiabatic sound speed and density) based on contemporary helioseismic 
measurements. They concluded that the mass losing Sun models were consistent with helioseismic 
measurements, with the 7% more massive sun case marginally more consistent than the standard 
solar model. While promising, the technique used by Sackmann & Boothroyd can neither support 
nor rule out a solar model with a mass loss of about 7% until improvements are made in helioseimic 
observational sensitivity and models of solar input physics. 
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The Sun currently loses a small amount of mass due to coronal mass ejections (~ 10 15 M Q yr 1 ] 
the solar wind (~ 2-3 x 1CT 14 M yr" 1 ) and radiation (~Tx 1CT 14 M Q yr" 1 ), the latter due to 



conversion of matter into energy through thermonuclear fusion at the core (jFeldman et al.lll977l ; 
Shulll982j ). Assuming that each of these rates have remained constant over time, the total amount of 
mass lost by the Sun over its 4.56 Ga history is only about 0.05% of the total solar mass. However, 
the physical processes that set the rate of mass loss due to solar wind and coronal mass ejections 
are poorly understood and it is not unreasonable to explore solar models in which the rates of mass 
loss were larger in the past by orders of magnitude. Astronomical observations indicate g enerally 
larger mass loss rates by stellar winds for younger solar analogues (jWood et al.ll2002j . 120051 ). There 
is some evidenc e in lunar regolith sa mples of a long term secular decrease in solar wind flux during 



the past 3 Gy (jKerridge et al 



1991 



There is also evidence in the meteoritic record of enhanced 
solar activity during the solar system's early history, though this has been interpreted as evidence 
of the Sun's very active but brief T Tauri phase, rather than an extended period of enhanced solar 



wind (jCaffee et al.lll987T ). 



In this paper we examine anew the hypothesis of a mass losing Sun. First we constrain the solar 
mass loss rate time history that would be required to solve the faint early Sun paradox by requiring 
that Earth's mean surface temperature remain above 273 K for all of its history (Section[2]). We then 
compare the resulting mass loss functions with estimates in the literature for the stellar wind mass 
loss rates of younger sun-like stars (Section [3]). The comparison is discouraging: the astronomical 
data does not support a solar mass loss history that would resolve the faint early Sun paradox. 
However, astronomical data cannot rule out the hypothesis either. Therefore, with a view towards 
identifying other tests of the hypothesis, we explore the effects of a mass losing Sun on the orbital 
dynamics of the solar system, including how a 7% more massive Sun would affect mean motion 
and secular resonances, the irregular satellites of the Jovian planets, and other small bodies in the 
solar system (Section H]). We summarize our results and conclusions in Section [5j 



2. Effect of Solar Mass on Climate History 



A very simple climate model can be used to explore the dependence on the mass of the Sun of 
the mean surface temperature of the Earth, while keeping all other factors constant. The steady 
state mean surface temperature o f the Earth (or any other terrestrial planet) is given by the 
following energy balance equation ()Pollacklll979l ): 



[l- A) SttR 2 = aeT^-KR 2 



(1) 



where A is the average planetary albedo, S is the solar flux at the top of the atmosphere, R is the 
planetary radius, e is the atmospheric IR emissivity, and a is the Stefan-Boltzmann con stant of 
5.67 x 10- 5 



a 



gm • s 



K . Typical values for the Earth are A = 0.34 and e = 0.6 (| Pollack 



19791 ; lHartmannlll994l ). For the purposes of this paper, we make the minimalist assumption that 
these parameters are time-independent and we take the current typical values to be constant over 
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all of geologic time. This is not necessarily the most accurate way to model the Earth's atmosphere, 
since changes in atmospheric greenhouse gas composition can vary the emissivity parameter e and 
changes to surface composition and cloud cover can change the planetary albedo parameter A; this 
parametrization may also be criticized as insufficient to represent the variety of greenhouse effects 
that may be possible. However our purpose is to identify solar mass histories that simply and 
completely resolve the faint early Sun paradox without the need to invoke other effects. Thus, for 
our purposes Equation CD provides an adequate representation, and can be rearranged to give the 
terrestrial surface temperature as a function of solar flux: 

nV4 



SQ.-A ) 



(2) 



Goughl (|198ll ) reports that, based on stellar nucleosynthesis and stellar evolution models, the 
rate of increase in luminosity of the Sun with time, assuming the solar mass is constant, can be 
represented as 

r o i- 1 

t/ t ) L , (3) 

where t = 4.56 x 10 9 yr is the current a ge of the Sun ( Dearbornl Il991h . and L Q ~ 3.9 x 
10 33 erg sec 



L(t) 



1+ | P 



= 4.56 x 10 9 yr is the current a ge of the Sun ( Dearbornl Il9£ 
1 is the current solar luminosity ( Shulll982l ). We will assume that the Earth's orbital 



eccentricity can be neglected and that the Earth's orbital radius is equal to its semi-major axis, a. 

A mass losing Sun will affect the Earth's climate in two ways. First, the main sequence stellar 
luminosity is quite sensitive to stellar mass, L oc M p , where P is in the range 3-5 depending upon 
source s of opacity (jib enl 1 196 71 ); for the early Sun, we adopt P = 4.75, following IWhitmire et al 



()1995l ). Second, if the solar mass loss is slow compared with the orbital motion of the planets, then 
it follows from the adiabatic invariance of the actions for a keplerian orbit that 



[M(t) + Mi] ai(t) = constant, 



(4) 



where Mj and a\ are the mass and orbital semimajor axis of a planet (or minor planet) in the solar 
system, and M(t) is the time- varying solar mass. In the solar system, Mj <C Mq, therefore a 
planet's semi-major axis will be smaller for a more massive Sun by a, oc M _1 . The solar radiation 
flux at Earth, S, is related to the Sun's luminosity L and Earth's semi-major axis a by S oc La~ 2 . 
Therefore, the time-variation of the solar flux at Earth is given by 

-1 / itrfj-W 6.75 



S(t) = So 



t/ t ) 



(M(t)Y 
V 



(5) 



where So is the current solar radiation flux at the top of Earth's atmosphere, So = 1.37 x 10 6 erg • 
cm -2 • s _1 . Using Eqnuation ([5]) in Equation ([2]), we obtain the time dependence of the terrestrial 



surface temperature: 



T.(t) 



So (1 - A) 



1/4 



Aea [1 + 1(1 -t/ t )] 



(M{t)\ 



1.69 



(6) 
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Note that Equation [6] is obtained from a straightforward energy balance model. We now use 
it to construct the "minimum mass-loss" history for the Sun that keeps the terrestrial surface 
temperature T s (t) above 273 K for all time, and is also consistent with the current solar mass loss 
rate of M Q ~ —2 x 10 -14 M yr _1 . We do this as follows. We assume that the current solar 
mass loss rate has remained constant for the past 2.3 Gy, the t ime span over which t he standard 
solar model yields terrestrial surface temperature above 273 K (jSagan Mullen! 1 1972l ). For times 
prior to 2.3 Ga, we calculate the solar mass that would be required to just keep the terrestrial 
surface temperature at 273 K. Therefore, we specify the terrestrial surface temperature at t = 0, 
T s (0) = 273 K. Equation © then yields the initial solar mass, M(0) = 1.026 M . Also from 
Equation ©, we obtain the solar mass as a function of time that keeps the terrestrial surface 
temperature constant, T s (t) = T s (0) = 273 K: 



M(t) = 0.973 M Q 



1+ | P 



t/ t ) 



0.148 



the corresponding mass loss rate is given by 



M(t) 



-1.26 x 10 



-ii 



n -0.852 



l + -(l-i/ t ) 



M y" 1 . 



(7) 



(8) 



In equations ([6]) and ([8]) we have adopted the numerical values of A and e typical of Earth mentioned 
above, and t = 4.56 x 10 9 yr. By demanding continuity of the solar mass time variation and 
consistency with the present solar mass loss rate, we obtain the piecewise "minimum solar mass 
loss" model: 

rO.974 M [l + |(l-t/t )]°- 15 for t< 2.39 Gy, 
\ M + M (t - t ) for t > 2.39 Gy. 

This function is shown by the short dashed curve in Figure [TJ 

Beyond the minimum solar mass loss history, we consider a solar mass loss rate that has 
decreased exponentially with time. We construct this with parameters constrained to keep the 
mean surface temperature of Earth above the freezing point of water for all of the planet's history. 



A functional form of this mass loss rate can be specified as M(t) 
implies solar mass as a function of time as follows: 



c + d(e' 



-at 



-a tr. 



M(t) = M + C(t - t ) + D(e 



-Qt 



-a t N 



which 



(10) 



If we denote the initial solar mass, M(0) = rut M , then the constants C and D are as follows: 

C = k-i«n P - tt 'e 
D 



M + aDe: 
(m/-l) M + M t 



1 _ e -a t Q 

where M = M( t ) ~ —2 x 10 -14 M y _1 is the present solar wind mass loss rate (at t = t ). 
The time constant, a, is found by considering the requirement that the mean surface temperature of 
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Earth remained above the freezing point of water during all of Earth's history: T (t < to) > 273 K. 
This yi elds a = 1.32 x 10 -9 yr _1 . The parameter mj is a free parameter, except that, as discussed 



above, iKastina (|1988l ) argues that the solar flux cannot have been more than 10% higher at any 



point in the past since Earth has not lost its oceans due to a runaway moist greenhouse atmosphere 



This p rovides an upper limit on the initial mass factor of mj = 1.07, and lSackmann Boothrovd 



(|2003l ) calculate that this value is consistent with helioseismology. With this value, we have C ~ 



-9.2 x 10~ n M y _1 and D ~ 0.070. The resulting solar mass function M(t) is shown in Figured! 

For the mass loss functions considered above, the corresponding mean surface temperature of 
Earth as a function of time is plotted in Figure [2j Note that although the cumulative mass loss of 
the Sun required to resolve the faint early sun paradox is modest, the timescale required for this 
mass loss is quite long, O(10 9 ) yr. 

We briefly consider Mars. With the standard solar model, in order to keep the global surface 
temperature on Mars above 273 K prior to 3.5 Ga, the effective atmospheric emissivity would need 
to be e = 0.29 (assuming its present albedo, A = 0.16), or e = 0.23 (assuming an Earth-like albedo, 
A = 0.34). For the mass losing solar models considered here, for the minimum required M that 
solves the faint early sun problem for Earth, in order to solve the problem for Mars the effective 
emissivity of Mars would need to be e = 0.33 (assuming its present albedo of A = 0.16) or e = 0.26 
(assuming an Earth-like albedo, A = 0.34). 



Kasting] (|199ll ) calculated what the required solar luminosities to keep Mars above the freezing 
point of water with a massive CO2 atmosphere and including the effects of CO2 condensation. 
He found that in order to keep early Mars warm enough for liquid water with a massive CO2 
atmosphere, the ratio of the solar radiation flux to the current solar radiation flux (S/Sq) had to 
be greater than 0.80-0.86. If the solar flux ratio was lower than this value, then the decrease in 
the convective lapse rate in the troposphere of Mars due to the condensation of CO2 clouds would 
offset the greenhouse warming enough that the surface could never be above the freezing point of 
water, no matter how high the CO2 pressure was. Figure [3] shows the solar radiation flux ratio as 
a function of time for the solar mass loss models considered here. Both the models considered here 
that solve the faint warm early Earth problem are also capable of satisfying Kasting's criteria for 
keeping early Mars warm. 



3. Stellar Winds of Sun-like Stars 

Since the sun is thought to be a typical G type dwarf star, it would be informative to measure 
the stellar wind mass outflows of other sun-like stars at various stages of their evolution on the 
main sequence. Stellar winds from G and K type dwarf stars are very difficult to observe directly 
owing to their low optical depth. Recently however, indirect measurements of the stellar winds 
of a small number of nearby sun-like stars have been made using the Hubble Space Telescope. 
These measurements exploit the charge exchange that occurs when the ionized stellar wind collides 
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with the neutral interst ellar medium, and which is detectable in a Hi Lya absorption feature (jZank 
19991 : IWood et al.ll2002l ). By measuring amount of Lya absorption in a star, and fitting the observed 
absorption feature to a model astrosphere interacting with its local interstellar medium, the mass 
loss rate of the star due to stellar winds can be estimated. The density and relative velocity of the 
local interstellar medium with respect to the star must also be known. 



Using this method, IWood et al.l ([20051 ) have estimated the stellar wind fluxes of several nearby 
solar-type main sequence stars. For each star, Wood et al. reported a mass loss rate due to stellar 
winds relative to the currently measured solar wind mass loss rate. They also reported the stellar 
x-ray luminosity, and the stellar surface area. They then us ed a power law relationship between 
the stellar x-ray flux and the age as reported by lAvresI (|1997l ). 



-1.74±0.34 



(11) 



We use Equation [HI together with the Sun's age as a reference point, to obtain an age estimate for 
each of the stars in Wood et al.'s sample. We also calculated the errors in the stellar ages using the 
error in the power law relationship between stellar x-ray flux and age. In addition, we also searched 
the literature for independent measurements of stellar ages, if available. The ages estimated from 
Equation [11] do not always agree with the ages found by other methods, and some are well outside 
the error bars of the x-ray flux ages. We consider here a sub-set of Wood et al's stars, i.e., only the 
eleven G and K stars for which mass loss rates were obtained. This data is listed in Table Q] and 
plotted in Figure 0[ Mass loss function that fit to the data of Wood et al. with a power law are 
also plotted in Figure HI along with the two solar mass loss history models described in Section [2j 

It is clear that the mass loss rate functions that we derived in the previous section to solve 
the faint early Sun paradox are inconsistent with the astronomical data on the mass loss rates of 
sun-like stars as found by Wood et al. The required solar mass loss rates at early times (0-2 Ga) are 
about 1-2 orders of magnitude higher than that for the observed solar analogues. The cumulative 
total stellar mass loss inferred from the Wood et al. data is < 0.003 Mq, in contrast with the 
required 0.03-0.07 M . The data also indicates a timescale of O(10 8 ) yr for the decrease of stellar 
wind flux at early times, whereas our hypothesis requires a slower decay, on timescale ~ 10 9 yr. To 
summarize, the faint young sun paradox requires a cumulative solar mass loss and early solar mass 
loss rates that are about 1-2 orders of magnitude and about one order of magnitude, respectively, 
higher than the astronomical data have revealed. The caveat is t hat th e stellar wind flux estimates 
of the younger solar analogs are model-dependent. IWood et al.1 (120021 ) estimate that there may be 
a factor of two uncertainty in their estimated mass loss rates due to the uncertainty in the stellar 
wind velocity. 



4. Dynamical Effects of Higher Mass Sun 



In the above analysis we have shown that in order for a mass losing Sun to solve the faint early 
Sun paradox, the Sun's mass had to be significantly higher for a period of time on the order of 2 Gy. 
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In this section, we explore what effects such a higher mass Sun might have had on the dynamics 
of the solar system, and whether some signatures of the early massive Sun may remain imprinted 
on the orbital dynamics of the present planetary system. A non-exclusive list of dynamical effects 
is discussed below. 



4.1. Planetary Orbits 



A slow solar mass loss conserves all the three actions of a keplerian orbit. Therefore, orbital 
eccentricities and inclinations of all planetary orbits are unchanged by a solar mass loss, and it 
was shown above that the semimajor axes of each of the planets are related to the time-varying 
mass of the sun through Equation HJ This means that, for an initial mass of the Sun greater than 
present by a factor mj, the net expansion of all planetary orbits has been a factor nif relative to 
the initial orbits, and the orbital periods T, L have increased by a factor m|. Thus, for rrif ~ 1.07, 
and considering the mass loss histories described in Section [21 the year would have been shorter by 
2-6% during the Archean epoch (Earth age 0.5-2.5 Gy, see Figure [1]). Indicators of the length of 
the year in the geological record (tidal rhythmites, banded iron form ations) provide poor precision 
for pre-Cambrian epochs (|Williamsll2000l : Eriksson &: Simpsonll2000i ). and are unlikely to provide a 
test of a shorter year . 

The slightly more compact ancient planetary system would have experienced relatively stronger 
orbital perturbations arising from the mutual gravitational interactions of the planets. The rates of 
secular precession of apsides and of orbit poles, which are proportional to Mj/Ti M Q , would have 
been faster by a factor rut. Although these precessional motions have b een linked to climate cycles 
on timescales of O(10 5 -10 6 ) years i n recent Ear th history (lHinnovll2000l ). the much lower resolution 
of the very ancient climate record (|Erwinll2006l ) makes it unlikely that it could be useful in testing 
these implications of the higher ancient solar mass. 



4.2. Minor Planet Resonance Locations 

Locations of mean motion resonances depend on the ratios of the semimajor axes between the 
planets and a test particle. It was previously shown that the semimajor axes of the planets are 
related to the time varying mass of the sun through Equation |4j For a Sun with a higher mass given 
by the factor m/, the change in the ratio of the semimajor axis of Jupiter to a massless particle 
can be found by: 

Qp/aj = mf (M + Mj) 

a' p /a'j ~ m f M e + Mj 1 ' 
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Where the primed values refer to the pre- mass loss semimajor axes. For a mass loss of 7%, this 
gives a change in the ratio of semimajor axes of 



a P /aj 
a'/a'j 



1.00006 



The assumption that cij cx Mg 1 is therefore a good one when considering the locations of mean 
motion resonances. It follows that the relative locations of mean motion resonances of the major 
planets remain nearly unchanged as the solar mass changes, however the strengths of resonant 
perturbations are not invariant: the fractional amplitude of resonant perturbations is proportional 
to ~ (Mi/ Mo) 1 / 2 , go that resonance widths (measured as a range in semimajor axis), Aa, scale 
as ~ rrij . The Kirkwood Gaps and the Hilda group of asteroids in the main asteroid belt may 
have been subjected to the effects of an expanding resonance width due to solar mass loss. For 
the Kirkwood Gaps, asteroids from near the edges of the resonances would be destablized, but 
this would not leave a trace of the original boundaries of the gaps in the present asteroid orbital 
distribution. For the surviving Hilda asteroids, the result would be a more compact final orbital 
element distribution compared to their initial orbits; the net changes estimated are small, not 
inconsistent with the observed population, but also consistent with other dynamical processes in 
the ear ly solar system (e.g ., resonance sweeping that accompanied the orbital migration of the 



planets iFranklin et al.l 12004 ) . 



One of the strongest resonant perturbations on asteroids is the vq secular resonance which 
defines the inner edge of the present-day asteroid belt. Considering the effects of the time varying 
solar mass, we find that the free precession rates of minor planets as well as the secular frequencies 
of the major planets are both proportional to mj, therefore the locations of secular resonances are 
invariant relative to the orbits of the major planets. Thus, we expect no signature of the solar mass 
loss in the inner boundary of the asteroid belt. 



4.3. Irregular satellites of the outer planets 

The irregular satellites of the outer giant planets are small objects having semimajor axes, 
eccentricities, and inclinations much higher than the regular satellites; the majority of them are on 
retrograde orbits. They are thought to be objects that formed in heliocentric orbits which were 
subsequently captured into stable orbits around the giant planets. 

The irregular satellites all orbit well within their planet's Hill sphere. In the restricted three 
body problem, the Hill sphere is defined by the radius 

rH=a2 U^J (13) 



where m c is the central mass, and rri2 is the large primary mass. If the solar mass were higher in 
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the past by a factor mf, then the Hill radius of each planet would have been 

-4/3 



m 



f 



(14) 



where r# is the present value. Thus, a mass losing Sun causes the Hill sphere of each planet 
to expand; this offers a novel way for the giant planets to capture satellites. This mech anism is 
similar to the pull-down capture mechanism proposed by iHeppenheimer fc Pored (119771 ). which 
invokes the expansion of the planet's Hill sphere during the gas accretion phase; however, the 
gas accretion phase is very short-lived, whereas the solar mass loss would occur over a much 
longer time scale. As the sun loses mass, its Hill's radius increases as shown in equation [TH 
Objects that are weakly bound to a planet, or are unbound but near the we ak stability boundary 
of th e planet may become permanently captured as the Sun's mass decreases (jAstakhov Sz Farrellv 
20041 ) . The efficiency of this mechanism is only weakly dependent on planet mass. We speculate 
that it may offer an explanation for the o bserved similarity in the i r regular satellite populations 
of Jupiter, Saturn, Uranus and Neptune (| Jewitt fc Sheppardl 120051 ; ISheppardl 120061 ) . The long 
timescale of this process is also an attractive means for evolv ing irregular satellites into secular 
resonances (cf. ISaha &: Tremaind Il993l : iNesvorny et al.1 l2003l ) . A detailed examination of these 
possibilities is beyond the scope of the present paper, but we hope to address them in future work. 

We note briefly that the growth of the giant planets as they accreted large amounts of gas 
during their formation also causes an e xpansion of the Hill sphere. This has been invoked in 
the "pull-down" mechanism proposed by IHeppenheimer &: Pored (|1977l ) for the capture of irregu- 
lar satellites. In this mechanism, the planetary Hill sphere expands by a much larger magnitude 
and over a much shorter period of time compared to the solar mass loss mechanism. Another 
mec hanism for changing pl a netary Hill spheres is the late orbital migration of the giant plan- 



ets (Hahn Malhotra 



1999 



Tsieanis 



irregular satellites (|Brunini Conicetill995l ) 



ct al 



2005) which could also have facilitated the capture of 



5. Conclusion 

We have calculated that the minimal cumulative mass loss of the Sun that would resolve the 
faint early Sun paradox is ~0.026 M . Models with a Sun which is up to ~7% more massive than 
present are also consistent with helioseimological constraints on solar interior evolution and with 
geological constraints. However, an important conclusion of our study is that the solar mass loss 
history that would resolve the faint early Sun paradox requires the Sun to remain moderately more 
massive than present for 1-2 Gy in its early history (Figure [4]) . 

Astronomical evidence (also shown in Figure U]) suggests that both the cumulative mass loss 
and the timescale for the mass loss (by means of stellar winds) of sun-like stars is significantly short 
of what is required to resolve the faint early sun paradox. The cumulative mass loss of solar age 
sun-like stars on t he main sequence, based on the power law mass loss rate time history given by 
Wood et al.1 (|2005l ). is <0.3%, and the mass loss rate declines rapidly in ~ 10 9 yr. 
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Of course, the astronomical data does not directly rule out the possibility that the Sun had 
a different time history of mass loss than that indicated by the compilation of stellar wind flux 
estimates of an ensemble sun-like stars of various ages. We therefore examined the effects of a solar 
mass loss history on the orbital dynamics in the planetary system, including possible signatures in 
the geological and climate record. Most effects are found to be too small to provide unambiguous 
tests of the early-massive-sun hypothesis. This hypothesis may offer a novel explanation for the 
capture of irregular satellites of the outer planets, including the capture into secular resonance for 
some Jovian irregular satellites. However, even if the current satellites of the giant planets preserved 
the signature of a more massive sun, other dynamical events in the early Solar System, such as the 
accretion and orbital migration of the giant planets may also provide viable explanations. With 
the recent accumulation of considerable data on the properties of irregular satellites, it may be 
possible to test the viability of these alternative mechanisms for the capture of irregular satellites. 
This requires more detailed theoretical development of the alternative scenarios for comparison 
with observations. 

We thank Alex Pavlov for helpful discussions on Mars climate issues. RM thanks Jack 
W. Szostak for inspiring her interest in this problem. We are grateful for research support from 
NASA's Origins of Solar Systems Research Program and from the NASA Astrobiology Institute. 
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Fig. 1. — Solar mass over time for the three different solar mass history models considered. 
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Fig. 2. — Estimated mean surface temperature of the Earth calculated using Equation [6] for the 
three different solar mass history models considered. 
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Fig. 3. — Solar radiation flux, S(t), for any planet relative to the current solar radiation flux, Sq. 
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Fig. 4. — Mass loss rates due to stellar wind of the G-K main sequence stars from lWood et al.l (120051 ) 
are plotted as black circles along with the three different solar mass history models considered. The 
vertical error bars are calculated assuming the measured stellar wind mass loss are known to within 
a factor of two. 
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Table 1. Estimated stellar wind mass loss rates of sun-lik e stars with measured astrospheres, 

from lWood et al.1 ifcoOfih . 



Name 


X-Ray Age a 


Independent Age 


Stellar wind M 




(xlO 9 yr) 


(xlO 9 yr) 


(xl0~ 14 M - yr" 1 ) 


Sun 




4.56±0.1 b 


2.5 ± 1.5 


a Cen 


3.4 ±0.2 


6.52 ± 0.3 C 


5.0 ±3 


e Eri 


0.72 ±0.3 




75 ±45 


61 Cyg A 


1.9 ±0.3 




1.3 ±0.8 


e Ind 


2.3 ±0.3 


1.3 ± 0.3 d 


1.3 ±0.8 


36 Oph 


0.88 ±0.3 




38 ±23 


A And 


0.36 ±0.2 




13 ±8 


70 Oph 


0.90 ±0.3 


0.50 ± 0.22° 


250 ± 150 


f Boo 


0.44 ± 0.2 


0.06 ± 0.03° 


13 ±8 


61 Vir 


6.5 ±0.6 


6.6 ± 3 d 


0.75 ± 0.45 


5 Eri 


2.0 ±0.3 




10.0 ±6 


DK Uma 


0.36 ±0.2 




0.38 ± .23 
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